Constraining the nuclear pairing gap with pairing vibrations 
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Pairing interactions with various density dependencies (surface/volume mixing) are constrained 
with the two-neutron separation energy in the Tin isotopic chain. The response associated with 
pairing vibrations in very neutron-rich nuclei is sensitive to the density dependence of the pairing 
interaction. Using the same pairing interaction in nuclear matter and in Tin nuclei, the range of 
densities where the LDA is valid in the pairing channel is also studied. 
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I. INTRODUCTION 

Studies on pairing effects in both nuclear matter and fi- 
nite nuclei have known intensified interests in the recent 
years There are two main approaches for pairing, 
depending whether the mean field is based on Gogny fi- 
nite range interaction or on Skyrme interaction. In the 
first approach, a similar functional is used in both the 
particle-hole channel and the pairing channel, although 
interactions are not exactly the same due to the density 
dependence of the pairing interaction. In the Skyrme ap- 
proach, the functionals are meant to be different in the 
two channels, as witnessed for instance by their density 
dependence. The use of a different interaction in the 
particle-hole channel and in the pairing channel has been 
justified a decade ago @; this is for instance the case 
of employing the Skyrme interaction in the particle-hole 
channel and a zero range density dependent interaction 
in the pairing channel. We shall focus on the Skyrme 
approach: in this case the pairing density functional is 
difficult to constrain and it has not been possible to de- 
rive an universal pairing interaction during past decades, 
using for instance the odd-even mass staggering on fi- 
nite nuclei (see e.g. [!, [3). This may indicate the need 
for another approach, using additional constrains: should 
the pairing density functional be extended, and are there 
additional relevant observablcs to constrain it ? 

Nuclear matter could help in constraining the pairing 
functional. This requires however to bridge nuclei and 
nuclear matter through LDA in the pairing channel: its 
condition of validity should be more systematically anal- 
ysed. It has been recently shown that the two paired neu- 
trons are spatially localised in low density medium which 
corresponds to the surface of the nucleus [j| . The same 
conclusion is drawn by analysing the di-neutron config- 
uration in the excited states @, Q , and also performing 
calculations in low density matter in [H, Q , mainly renew- 
ing the possibility to link in some cases the nuclear matter 
and nuclei in the pairing channel through the LDA. 

Concomitantly the pairing functional has been ex- 
tended in order to study the condensation of the Cooper 
pairs (BEC-BCS crossover) in both symmetric and neu- 
tron matter. In nuclear matter the medium polarization 
increases the pairing gap at low densities in symmetric 
matter, whereas it reduces the gap in neutron matter, in- 



dicating an isospin dependence of the pairing functional 
[lOj , The application to finite nuclei of extended pairing 
density functional have shown the relevance of the LDA 
in the pairing channel (Tlj . 

The pairing functional studies may thus enter in a new 
era, renewing the method to design the pairing interac- 
tion: i) using an isospin dependence of the pairing inter- 
action ii) using eventually the nuclear matter as an addi- 
tional constrain for the pairing interaction iii) looking for 
additional observables in nuclei than the odd-even mass 
staggering to constrain the pairing interaction. Point i) 
has been investigated in [l(J, [ll[ . Point ii) requires the 
validity of the LDA in the pairing channel and shall be 
addressed in the present work. 

In the case of point iii) an interesting observable is pair- 
ing vibrations, measured through two-particle transfer. 
It is well known that the transfer cross section crucially 
depends on the pairing interaction at work in the trans- 
ferred paired [l2|, EH- However in the 70-80's the form 
factor of the transition has never been calculated fully 
microscopically. The first microscopic calculations has 
been performed only recently (l4| . allowing for a strong 
link between the pairing interaction and pairing vibra- 
tions. Several calculations followed [1, [HI, showing the 
renewed interest for such studies. 



It is therefore meaningful to use pairing vibrations as 
a complementary observable to the masses, in order to 
constrain the pairing interaction, and study the implica- 
tions to the nuclear matter. More generally, the isospin 
extension of the pairing interaction should be accompa- 
nied with additional constraints. One purpose of this 
work is to evaluate if pairing vibrations could play this 
role (Section III). 

The method is to analyse the sensitivity of pairing vi- 
brations to various pairing interactions which provide the 
same two-neutron separation energy in Tin isotopes, and 
evaluate the consequences on the pairing gap in symmet- 
ric and neutron matter. On this purpose it is necessary 
to determine the range of density where the LDA is valid 
in the pairing channel (Section II). 
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TABLE I: Values of r\ and Vo of the pairing interaction. 



0,16E 



n 


Vb (MeV fm- 1 ) 


0.35 


-285 


0.65 


-390 


1 


-670 



II. VALIDITY OF THE LDA IN THE PAIRING 
CHANNEL 



After many years of study, there is still no unambigu- 
ous universal pairing functional ranging on the whole 
nuclear chart, and current efforts are aiming in that di- 
rection. The problem may be due to the method used 
to constrain it, namely comparing the pairing gap with 
odd-even mass differences, or evaluating the separation 
energies along a given isotopic chain. It therefore may be 
useful to consider a more general context: the evaluation 
of several pairing interactions constrained by odd-even 
mass difference, on nuclear matter on one side, and on 
additional obscrvablcs on the other side, should shed a 
renewed light on the problem. To achieve this goal it is 
first necessary to determine the range of validity of the 
LDA in the pairing channel. 



Method to determine the functional 
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FIG. 1: Upper panel: Matter densities calculated with the 
HFB model for 124 Sn and 136 Sn. The vertical lines indicate 
the radius corresponding the density at which all the pairing 
interactions converge in uniform matter (see text). Lower 
panel: pairing field of 124 Sn calculated with a surface 77 = 1 
(top) and a mixed r\ = 0.35 (bottom) pairing interaction. 



As an illustration to visualize the features of the cal- 
culated pairing effects, we display in Fig. [T]the neutron 
pairing field for 124 Sn corresponding to the surface rj = 1 
and the mixed 77 = 0.35 interactions. 



The method is the following: we first consider surface 
and various mixed paring interactions. The parameters 
are determined so as to describe the two neutron sepa- 
ration energy. Then pairing vibrations are used in order 
to disentangle between the various pairing interactions 
(Section III). We choose 124 Sn and 136 Sn nuclei: these 
are spherical nuclei where pairing vibrations arc likely 
to occur [l3j]. One is stable and the second has a large 
neutron excess. 

The microscopic calculations for the ground state are 
based on the Hartrec-Fock-Bogoliubov (HFB) model. 
The Skyrmc interaction SLy4 [l6[ is chosen for the 
particle-hole channel of the HFB equations. The adopted 
pairing interaction is the usual zero-range density- 
dependent interaction 



Vp a i r — Vq 



p( r ) 

PO 



5 (rr - r 2 ) 



(1) 



where 77 provides the surface/volume character of the 
interaction. We set a = 1 and po — 0.16 fm -3 . The nu- 
merical cutoff for the microscopic calculations is given by 
Emax = 60 MeV (in quasiparticle energies) and j ma x— 
15/2. For each value of 77, Vb is chosen to fit the known 
experimental two-neutron separation energies for Sn iso- 
topes. Surface and mixed interactions have been consid- 
ered in this work and the used values of (77, Vb) are listed 
in Table I. 



B. Pairing gap in uniform matter 

The relation between the pairing gap in uniform matter 
at a given density and the pairing field at a given radius 
in nuclei has been explored in Ref . [ll[ . It has been found 
that in the case of mixed interactions, the LDA is in good 
agreement with the full microscopic HFB calculation (dif- 
ferences less than 15% on the pairing field). This might 
be related to the extension of the Cooper pair which is 
getting smaller at the surface of nuclei (about 2 fm) com- 
pared to that in the interior (about 5-6 fm) [f| . Close to 
the surface, pairing properties shall not be very differ- 
ent from that of a uniform piece of matter at the same 
density. It is then interesting to explore the low density 
properties of the different pairing interactions listed in 
Table I. 

Fig. [2] displays the pairing gap in uniform matter for 
various pairing interactions. It is observed that the differ- 
ent interactions leads to very different pairing gap at low 
density while around saturation density, there is a den- 
sity (p=0.11 fm -3 ) at which the pairing gap and pairing 
strength coincide for the three pairing interactions. 

From Fig. [21 two conclusions can be drawn: i) the two- 
neutron separation energy used to adjust the parameters 
of the pairing interaction is an observable which provides 
a strong constrain on the pairing gap localized at the sur- 
face of the nuclei: the pairing gap in nuclear matter has 
been constrained for p=0.11 fm -3 , which corresponds to 
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FIG. 2: Pairing gap versus the density for uniform matter for 
different pairing interactions. 

R ~ 5 fm in Tin nuclei, as showed by the 124 Sn densities 
displayed on Fig. [T] ii) to better constrain the value of 
the parameter 77, one shall find another observable sensi- 
tive to the pairing strength at low density (large radius, 
experimentally easier to probe). Indeed, in the very ex- 
ternal part of the nuclei the pairing strength is very dif- 
ferent from one interaction to another. The pure surface 
pairing interaction predict a pairing gap as high as 8 MeV 
at low density while the various mixed pairing interaction 
are grouped below 3 MeV (see Fig. [2|). 

Therefore, one might expect that properties of collec- 
tive modes sensitive to the external part of the nuclei 
could be changed by the properties of the pairing inter- 
action at low density. Pair transferred reaction mecha- 
nisms like (p,t) or (a, 6 He) which is very surface peaked 
shall also help in extracting the value of the pairing gap 
in the external part of the nuclei or equivalcntly at low 
density. 

C. Validity of the LDA 

One aim of this work is to investigate the mapping 
of the pairing gap between nuclei and nuclear matter, 
implying the following relation: 

Aoo(p) = A HFB (p(r) = p) (2) 

Eq. (J5J) corresponds to the LDA in the pairing channel. 
However, it should be noted that in practice, the LDA 
might still be valid even if Eq. ([2|) is not verified for very 
low densities (p < 0.1 fm -3 ): the major role is played by 
densities around the saturation density, where nucleus 
properties are the most important. 

Fig. [3] displays the gap calculated with the HFB ap- 
proach as a function of the nuclei density in 124 Sn. Nu- 
clear matter gap is surimposed, where the neutron-proton 



FIG. 3: Pairing gap versus the density for nuclear matter and 
for 124 Sn (see text) 

asymmetry are taken as the one of 124 Sn. One can see 
that the LDA is not valid for a pure surface pairing, but 
becomes valid in the 77=0.35 mixed case, for densities 
lower than 0.1 fm -3 . Therefore, in this case, the low den- 
sity part of the nuclear matter gap could be constrained 
by nuclei calculations of the gap: a valid case of the LDA 
in the pairing channel is showed here quantitatively for 
the first time. This range of density corresponds to a 
radius larger than about 5.1 fm in 124 Sn and 5.2 fm in 
136 Sn localised at the surface of the nuclei, as can be seen 
on Fig. Q] In the cases of the other pairing interactions, 
the LDA underestimates the pairing gap at low density 
and overestimate it at high density. However, it should 
be noted that the relative difference of the gaps between 
the various pairing interactions are respected. 

As noted above, to assess in practice the validity of 
the LDA, the density range around the saturation den- 
sity should be investigated, where nuclei bulk properties 
are at work. Therefore a large discrepancy at very low 
density (p < 0.1 fm -3 ), as observed in the case of the sur- 
face pairing interaction in Fig . [31 does not preclude on 
the validity of the LDA. More generally it can be stated 
that the less important the surface contribution to the 
pairing interaction, the better the LDA. 



III. PAIRING VIBRATIONS 

As stated above, it may be useful to consider an addi- 
tional observable than the separation energy, in order to 
constrain the pairing interaction, namely its density de- 
pendence. There are only few observables which could be 
relevant to constrain pairing effects. It has been shown 
that the first 2 + state in nuclei is sensitive to the pairing 
interaction [l7| : both its position and strength depend 
of the pairing interaction. However this is mainly related 
to the pairing gap value, which is the same observable 
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extracted from odd-even mass excess. It should be noted 
than none of these two obscrvablcs (the first 2 + state and 
the odd-even mass staggering) can be directly linked to 
predictions. On one side there is the difficulty to mod- 
elize excited states. On the other side, the difficulty is to 
describe odd nuclei. 

Pairing vibrations may be more adequate observable. 
They can be probed for instance with two neutron trans- 
fer in nuclei close to shell closure. We refer to [H, [H| 
for details on pairing vibrations. Basically, these modes 
corresponds to the (collective) filling of subshells, in tran- 
sition from an A to A+2 nuclei. 

With pairing vibrations, pairing effects are probed by 
3 ways. The first one is the magnitude of the pairing 
gap A (average of the pairing field): a large pairing gap 
implies strength at larger energies, following the formula 
E 2 ~(e — A) 2 + A 2 . This component is also present in 
the first 2 + state in the ph response as well as in the 
odd-even mass staggering. But in the case of the pairing 
vibrations, there are two additional contributions: first, 
the transition densities generating the strength are the 
pairing one, which means that the unperturbed response 
as well as the perturbed response are sensitive to the 
impact of the pairing on the wave functions. Finally, the 
residual interaction, generating the QRPA response, is 
also sensitive to pairing. 

Pairing vibrations are therefore expected to be very 
sensitive to the pairing interaction. On the other hand, 
it may also be difficult to disentangle between the three 
above mentioned effects. However, the first one can be 
evaluated using the energies of the unperturbed response, 
the second one by studying the pairing transition densi- 
ties, and the last one by comparing the unperturbed and 
the QRPA responses. It should be noted that a related 
study will also be performed in (l8| . 

A. Method: QRPA in the pp channel 

The method is described in [13, [13 ■ Namely the QRPA 
equations are solved in coordinate space, using the green 
function formalism. The variation of the generalized den- 
sity TV is expressed in term of 3 quantities, namely p' , n' 
and k', which are written as a column vector: 




where p[j = {0\cjCi\'J is the variation of the parti- 
cle density, k'^ = (0|cj-c,|') and k'^ = /o|c]cj|'^ are the 

fluctuations of the pairing tensor associated to the pair- 
ing vibrations and |'} denotes the change of the ground 
state wavefunction |0 > due to the external field. In con- 
trast with the RPA where one needs to know only the 
change of the ph density (//), the variation of the three 



quantities (J3)) have to be calculated in the QRPA. In the 
three dimensional space introduced in Eq. ([3]), the first 
dimension represents the particle-hole (ph) subspace, the 
second the particle-particle (pp) one, and the third the 
hole-hole (hh) one. The response matrix has therefore 9 
coupled elements in QRPA, compared to one in the RPA 
formalism. 

The variation of the HFB Hamiltonian is given by: 

H' = Vp', (4) 

where V is the matrix of the residual interaction ex- 
pressed in terms of the second derivatives of the HFB 
energy functional, namely: 

V^(ra,rV)= a, /3 = 1, 2, 3. (5) 

In the above equation the notation a means that when- 
ever a is 2 or 3 then a is 3 or 2. 

The QRPA Green's function G can be used for cal- 
culating the strength function associated with the two- 
particle transfer from the ground state of a nucleus with 
A nucleons to the excited states of a nucleus with A+2 
nucleons. This strength function is : 

S(uj) = --7m J F*(r)G 22 (r,r';Lj)F(r')dr dr' (6) 

where G denotes the (pp,pp) component of the Green's 
function and F in the external perturbating field associ- 
ated with the addition of two particles. 

In the QRPA calculations the full HFB quasiparticlc 
spectrum up to 60 MeV is included. These states are 
used to construct the unperturbed Green's function Go- 
The residual interaction is derived from the two-body 
force used in HFB according to Eq. ([5]). The con- 
tribution given by the velocity-dependent terms of the 
Skyrme force to the residual interaction is calculated in 
the Landau-Migdal approximation, which is shown to be 
accurate [l9j . The strength function for the two-neutron 
transfer is calculated using Eq. ([6]). The unperturbed 
Green's function is calculated with an averaging interval 
equal to 0.15 MeV. All details can be found in Ref. [iH ]. 

The response function is calculated for the pp channel. 
All the calculations are performed in a box of size 22.5 
fin. It should be noted that exact continuum treatment 
is much heavier, especially for nuclei such as Sn isotopes. 
Moreover the aim is not to study the impact of the con- 
tinuum treatment (see [l4| for such a study). Finally the 
Sn isotopes under study are far from the drip line, and 
continuum effects are expected to play a negligible role. 

B. Unperturbed response results 

The HFB solutions are used in the quasiparticlc 
random-phase approximation (QRPA) scheme to analyze 
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sclf-consistcntly the excitations modes associated to the 
pair transfer reactions. Since we study here two-neutron 
transfers, we focus on the neutron HFB quasiparticlc 
states that are used to construct the elementary config- 
urations of the excited modes. We work with positive- 
energy quasiparticle states. Once calculated the quasi- 
particle spectrum, it is possible to deduce some proper- 
ties of the unperturbed response function. 

The quasiparticle states with energy less than 6 MeV 
and an occupation probability < 80 % are presented in 
Tables II and III for 124 Sn and 136 Sn, respectively. Let 
us discuss the two cases r\ = 0.35 and r] = 1 (for T] = 0.65, 
results are similar to those obtained with rj = 0.35). For 
124 Sn, in the case of a mixed pairing interaction, r/ = 0.35, 
all the quasiparticle states with energy lower than 5 MeV 
are totally occupied with the exception of a /111/2 state 
at 1.5 MeV which is 42% occupied. This is the only low- 
energy state that can contribute to some extent to the ex- 
citation mode. The states that are completely empty and 
can thus contribute more to the excitation are located at 
higher energies. The first is an f-jj 2 state at 5.8 MeV. 
The others have larger energies (at least 1 MeV more). 
One can thus expect that the unperturbed response pro- 
file starts with a peak at twice 5.8 MeV, i.e. at ~ 11.6 
MeV (with some small contribution at 3 MeV). In the 
case of a surface interaction, rj = 1, again, all the states 
between and 5 MeV are occupied with the exception of 
a hn/2 state at 2.2 MeV (42% of occupation). This time 
there are several unoccupied states just above 5 MeV, the 
lowest energy being at 5.4 MeV (7)3/2 state). Hence, the 
unperturbed response is expected to have some structure 
starting from ~ 10.8 MeV with a small contribution at 
- 4.4 MeV. 

For the nucleus 136 Sn the situation is different: there 
are several low-lying unoccupied states. For rj = 0.35 the 
lowest energy for a completely unoccupied state is 1.9 
MeV (P3/2 state). At 0.8 MeV one also finds a fj/ 2 state 
with 45% of occupation. In the case 77 = 1 the lowest 
energy for a totally unoccupied state is 1.7 MeV (P3/2 
state) and a / 7 / 2 state is found at 1.6 MeV with 32% 
of occupation. The unperturbed response is expected to 
start at ~ 3.8 and 3.2 MeV for r) = 0.35 and 1, respec- 
tively. In the former case a small contribution at ~ 1.6 
MeV is also expected. 

In order to disentangle between the various pairing ef- 
fects, the unperturbed response in the two neutrons ad- 
dition mode is first shown on Fig. [4] for 124 Sn. The un- 
perturbed response is built on the HFB single quasipar- 
ticle (QP) spectrum, for the three pairing interactions. 
It should be noted that the spectrum is showed above 10 
MeV, because there is only the h n / 2 subshcll which can 
welcome two neutrons to make a low energy state: all 
the other configurations belong to the next major shell 
(see Table II) , explaining this high energy feature of the 
spectrum, as stated above. For all the mixed pairing in- 
teraction, the unperturbed spectrum is similar, showing 
that both the single quasiparticle energy and wave func- 
tions are close to each other in that case. However, in the 



TABLE II: Neutron quasiparticle states with E < 6 MeV and 
occupation less than 80%. The nucleus is 124 Sn. 
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State 


E (MeV) 


occ 


0.35 


hll/2 


1.5 


0.42 




17/2 


5.8 


0.01 


0.65 


hll/2 


1.7 


0.42 




f7/2 


5.7 


0.01 


1 


hll/2 


2.2 


0.42 




P 3/2 


5.4 


0.003 




f7/2 


5.5 


0.02 




pl/2 


5.6 


0.002 




sl/2 


5.7 


0.002 




TABLE III: Same 


as in Table II but for 136 Sn. 


V 


idle 


E (MeV) 


occ 


U.OO 


f7 /9 

li j A 


0.8 


0.45 




t-»^ /9 
P6/Z 


1.9 


0.01 




pi/2 


2.4 


0.006 




Pi /9 


2.9 


0.01 




sl/2 


3.3 


0.0005 






4.0 


0.0002 




d3/2 


4.1 


0.0005 




rfQ /9 


5.6 


0.0001 




<r7/9 


5.6 


0.0001 


U.DO 


f7 /9 
It / z 


0.9 


0.43 




T~»9 /9 


1.9 


0.02 




pl/2 


2.4 


0.008 




Pi /9 


2.9 


0.02 




sl/2 


3.2 


0.0004 




d5/2 


3.9 


0.0001 




H 1 } /9 

OO / £* 


3.9 


0.0003 




o-Q/9 


5.4 


0.0001 




p-7/9 


5.5 


0.00003 


1 
1 


f7 /9 
II / Z 


1.6 


0.32 




t-. Q /O 
pj/2 


1.7 


0.02 




pl/2 


1.9 


0.01 




sl/2 


1.9 


0.0004 




d5/2 


2.6 


0.0003 




d3/2 


2.6 


0.0002 




f5/2 


3.0 


0.01 




g9/2 


4.1 


0.0002 




g7/2 


4.1 


0.0001 



case of the pure surface pairing, the spectrum is changed. 
The energies are shifted to lower values, and the overall 
strength is increased. The lower energy shift can be un- 
derstood by more single QP states located at low energy. 
This can be explained by a lower pairing gap and a dif- 
ferent energy spectrum found in the HFB self-consistent 
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FIG. 4: Unperturbed response function for 124 Sn in the two 
neutrons + addition mode. The pure surface mode is in solid 
line, the -q=0.65 mode is in dotted line, and the ?7=0.35 mode 
in dashed-dotted lines 



FIG. 5: Unperturbed response function for Sn in the two 
neutrons + addition mode. The pure surface mode is in solid 
line, the -q—0.65 mode is in dotted line, and the »7=0.35 mode 
in dashed-dotted lines. 



procedure. The larger magnitude comes from the wave 
functions, and will be studied with the QPRA response. 
It can be already stated that the QRPA response will also 
have more strength at lower energy, due to this peculiar 
feature of the unperturbed spectrum for the pure surface 
pairing force. 

Fig. [5] shows the unperturbed response for the two 
neutron addition mode in 136 Sn. In this case, at the 
beginning of an open neutron shell several low energies 
configurations can welcome the two neutrons (see Tabic 
III). As in the case of 124 Sn, the response exhibits larger 
strength at low energy in the specific case of the pure 
surface pairing interaction, compared to others pairing 
interaction. This is related to the pairing field profile as 
shown on Fig. [T] It should be noted that in order to 
clearly sec the effect due to the surface pairing, not only 
the first + state, but also the energy area of a few MeV 
above should be explored since the results are different 
from to 4 MeV on Fig. [5] 

C. Perturbed response results 

Fig. [6]shows the QRPA response for 124 Sn, with a pure 
surface and the two mixed interactions. As expected the 
residual interaction plays a similar role in all the cases, 
gathering strength and shifting it to lower energy. In the 
case of 124 Sn, a peak around 9 MeV is the strongest for 
the surface pairing interaction, to be compared with the 
one around 10 MeV for the other interactions. Hence it 



is expected that the pairing vibration transition strength 
should be larger in the case of a pure surface force. How- 
ever it is known that it is difficult to describe accurately 
the magnitude of these transitions, especially for abso- 
lute cross section calculations [2(J: one-step or sequen- 
tial two-step process, triton wave function, zero-range or 
finite-range DWBA have to be considered. It is therefore 
necessary to rely on the angular distribution, calculated 
from the form factor, related itself to the pairing transi- 
tion density [l3[ . 

The pairing transition density is defined as: 

K v (r,a) = (0\c(r,a)c(r,a)\u) (7) 

where (r,a)= —2<tc> (r, —a) is its time reversed coun- 
terpart. 

It allows to calculate the form factor in the zero-range 
DWBA approximation. The pairing transition densities 
of Fig. [7] show, in the case of 124 Sn, a difference, go- 
ing from surface to other modes: the transition density 
decreases at the surface. However the difference is not 
dramatic and may be overruled by the experimental un- 
certainties. The larger strength of the 9 MeV peak in the 
pure surface pairing interaction is due to a larger density 
at the surface. 

For the 136 Sn neutron-rich nucleus, the low energy 
spectrum displayed on Fig. [5] is dramatically changed 
from using surface to other interactions, on a several MeV 
area. A three peaks structure appears in the surface case, 
compared to the 2 peak structure of the other cases. The 
integrated strength is also larger in the surface case. 
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E* (MeV) E* (MeV) 



FIG. 6: QRPA response function for Sn in the two neutrons 
+ addition mode. The pure surface mode is in solid line, 
the 77=0. 65 mode is in dotted line, and the 77=0.35 mode in 
dashed-dotted lines. 



FIG. 8: QRPA response function for Sn in the two neutrons 
+ addition mode. The pure surface mode is in solid line, 
the 77=0. 65 mode is in dotted line, and the 77=0.35 mode in 
dashed-dotted lines. 



0.03 - 




r(fm) 



Fig. [9] and \TU\ show the corresponding transition den- 
sities. They exhibit very different shapes, comparing re- 
sults with the pure surface pairing interaction and the 
mixed pairing interaction. Hence 136 Sn is a good test 
case to probe the pairing interaction through pairing vi- 
brations. For instance in the case of the most intense 
peak, the central part is dominant in the transition den- 
sity for the mixed case, whereas the surface part of the 
transition density dominates in the pure surface interac- 
tion. Hence a measurement of the angular distributions 
associated with the pairing vibration strength in very 
neutron rich- nuclei such as 136 Sn seems more decisive to 
disentangle between the pairing interactions than with 
136 Sn. This may be due to the larger neutron skin in 
136 Sn, consisting of low density neutron-rich matter. 

It has been shown in a previous work how the pairing 
transition densities allow to calculate the two neutron 
form factor in order to predict angular distributions for 
pairing vibrations [3] . Work along these lines should be 
undertaken in order to bring additional constrains on the 
pairing interaction. 



FIG. 7: Neutron transition density in the two neutrons addi- 
tion mode for 124 Sn for the first peak located at 9-10 MeV. 
The pure surface mode is in solid line, the 77=0.65 mode is in 
dotted line, and the 77=0.35 mode in dashed-dotted lines. 



IV. CONCLUSIONS 

The impact of various pairing interactions on pair- 
ing vibrations predictions has been analysed for the first 
time using a HFB+QRPA approach. They should pro- 
vide a good sensitivity from a pure surface interaction 
compared to mixed interactions, especially in the case 
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of very neutron-rich nuclei such as 136 Sn. Moreover nu- 
clear matter gap calculations show that the LDA in the 
pairing channel is valid in the surface of the nuclei for 
the 77 = 0.35 mixed surface/volume interaction, but not 
rigorously valid in other cases. In the case of exotic nu- 
clei, pairing vibrations are found more sensitive to the 
surface/volume type of the pairing interaction, than in 
the case of stable nuclei. This may be due to the larger 
extension of the neutron density in very neutron-rich nu- 
clei. 



The same study using an isospin dependent pairing in- 
teraction will be undertaken. The hope is to come one 
step closer to a more global pairing interaction, using 
odd even mass staggering, pairing vibrations, and nu- 
clear matter as constraints. Experimentally, the pairing 
transition densities can be tested through the form factor 
used to calculate the two neutrons transfer cross section. 
This implies to use a adequate reaction model. Work 
along these lines will be undertaken in an near future. 



FIG. 9: Neutron transition density in the two neutrons addi- 
tion mode for 136 Sn for the first two peaks of the strength, in 
the case of the mixed 77=0.65 interaction 
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FIG. 10: Neutron transition density in the two neutrons ad- 
dition mode for 136 Sn for the first 3 peaks of the strength, in 
the case of the pure surface interaction 



